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Abstract 

The sinh-Gordon model on a half-line with integrable boundary conditions is 
considered in low order perturbation theory developed in affine Toda field theory. 
The quantum corrections to the classical reflection factor of the model are studied 
up to the second order in the difference of the two boundary parameters and to 
one loop order in the bulk coupling. It is noticed that the general form of the 
second order quantum corrections are consistent with Ghoshal's formula. 
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1 Introduction 



Affine To da field theory [[|, is one of the most successful examples of the two- 
dimensional quantum field theory since it can be solved exactly. This theory possesses 
remarkable properties including classical and quantum integrability. An interesting re- 
view of the recent developments in Affine Toda field theory is presented in Ref. 0. The 
classical affine Toda field theories remain integrable in the presence of certain bound- 
ary conditions restricting them to a half-line, or to an interval [f|-||. Indeed, Corrigan 
et.al have classified the boundary conditions which preserve classical integrabil- 
ity. However, quantum integrability are hardly explored in the presence of a boundary 
although there has been progress for models based on class of theories P, |IT)|-|i"3]] . 
The simplest of these (n=l) is the sinh-Gordon model. In fact, this model is the only 
example in the ade series of affine Toda field theory which allows continuous boundary 
parameters. Although the sinh-Gordon model has been studied much more than other 
models in the integrable quantum field theories there still remains much to be studied 
in connection with the quantum integrability of the model with a boundary. In par- 
ticular, one of the interesting problems is to see how the boundary parameters could 
relate to the quantum reflection factor of the theory. 

In 1993 Ghoshal and Zamolodchikov |TJj studied the behaviour of the sine-Gordon 
model restricted to a half-line. They found the soliton reflection factors in the model 
and discovered a two parameter family of boundary conditions which preserved inte- 
grability. Then, Ghoshal computed the reflection factors of the soliton-antisoliton 
bound states (the breathers) of the model. The pioneering ideas of Ghoshal and 
Zamolodchikov motivated many researchers to investigate integrable quantum field the- 
ory with a boundary. Corrigan [TE\ was the first to discover that the lightest breather 
reflection factor of the sine-Gordon model is identical to the reflection factor of the 
sinh-Gordon model followed by an analytic continuation in the coupling constant. 

Recently, Corrigan and Delius |T7j calculated the bound-state spectrum of the sinh- 
Gordon model on a half-line in two different ways, firstly by using a boundary bootstrap 
principal, and secondly by quantizing the classical boundary breather states using a 



WKB technique (for further discussion see Ref. [fL8|). Comparing these two calculations 
provides a relationship between the boundary parameters, the bulk coupling constant 
and the parameters in Ghoshal's formula. They performed the calculations in the 
special case where the bulk Z 2 symmetry is conserved at the boundary, by requiring 
the two boundary parameters to be equal. Then, Corrigan and Taormina [|19| studied 
and extended the idea of |T7| in the intricate case where the two boundary parameters 
are different and the bulk symmetry is violated. Their calculations lead to a determi- 
nation of the relationship between the Lagrangian boundary parameters and the two 
parameters appearing in the reflection factor of the sinh-Gordon particle. Moreover, 
their results clarify the weak-strong coupling duality of the sinh-Gordon model with 
integrable boundary conditions. 

The one loop quantum corrections to the classical reflection factor of the sinh- 
Gordon model have been found |I6] when the boundary parameters are equal. If the 
two boundary parameters are different and the bulk Z 2 symmetry of the model is 
broken then, the lowest energy static background configuration will not be a trivial 
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solution {(f) = 0). In this intricate case, the perturbation theory must be developed 
within a complicated background. Meanwhile, the corresponding calculations involve 
a significantly more difficult propagator as well as intricate coupling constants. In a 
recent paper |^D| the quantum reflection factor has been found in one loop order up to 
the first order in the difference of the boundary parameters. The calculations provide a 
further verification of Ghoshal's formula up to the first order. The quantum reflection 
factor up to the second order in the difference of the boundary data is treated in this 
article. The result could test more deeply the expression for the sinh-Gordon model 
reflection factor given in Ref. [TE . 



2 The sinh-Gordon model on the half-line 

The sinh-Gordon theory describes a real massive scalar quantum field theory in 1+1 
dimensions with exponential self-interaction. The Lagrangian density of the theory in 
the presence of a boundary is 

£ = 6{-x)C - 6{x)B (2.1) 
where the bulk Lagrangian density of the model is 

C = ~d^d»<i>-V(<f>) 

i 2m 2 
= cosh(V2/ty), (2.2) 



and the boundary potential B has the following generic form WA 



B=j 2 (a Q e-^ + a l e^). (2.3) 

Note, in the above expressions m and (3 are a mass scale and a coupling constant of 
the theory and we have used normalizations customary in affine Toda field theory. 
The two real coefficients <7o and o\ are essentially free which indicate |5],|7j the extra 
parameters permitted at the boundary x = 0. Meanwhile, the boundary potential is 
required to satisfy the following equation as a consequence of maintaining integrability 
on the half-line: 

i = -f ™ 

or 

^ = ( aie WV2 _ a e -WV2\ at x = o (2 5) 

ox p \ J 

where we have rescaled the mass scale to unity in the integrable boundary condition. 
It is also convenient to use an alternative expression for the boundary parameters, i.e. 

Oi = COSaj7T. 

The sinh-Gordon model is integrable classically which implies there are infinitely 
many mutually commuting, independent conserved charges Q± s , where s is an arbitrary 
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odd integer. On the other hand, the model is integrable in the context of quantum field 
theory which means the S-matrix describing the n-particles scattering factorises into 
a product of two-particles scattering amplitudes. The S-matrix describing the elastic 
scattering of two sinh-Gordon particles with relative rapidity 9 is conjectured to be 
given by H|l],|2 

S{B) = -(BW^Bj <2 ' 6) 
where we have used the hyperbolic building block notation 

sinh(#/2 + is) 
x) = _ 4 ; (o 7) 

[) mnh(e/2-ifY 

and the quantity B is related to the coupling constant j3 by B = ■ 

For the boundary sinh-Gordon model, the most important physical quantity is the 
boundary S-matrix or the reflection factor which describes the reflection of a single 
particle from the boundary. Firstly, Ghoshal and Zamolodchikov |TJJ calculated the 
soliton reflection factors for the sine-Gordon model by solving the boundary Yang- 
Baxter equation and using general constraints implementing unitarity and a form of 
crossing symmetry. Then, Ghoshal |T5|] calculated the soliton-antisoliton bound states 



reflection factor by using the boundary bootstrap equation. Finally, following the idea 
discovered by Corrigan f[6l , the quantum reflection factor for the lightest breather in 



the sine-Gordon model is supposed to be identical with the quantum reflection factor 
of the sinh-Gordon particle after doing analytic continuation in the coupling constant 
to derive 



R (9) = (l)(2-B/2)(l + B/2) 

q{ ' (1 - E(a , <r u /3))(1 + E(a , a u - F(a , a 1} /3))(1 + F(a , a u (3)) ' 

(2.8) 

When the bulk reflection symmetry is preserved which means do = <J\ , then one of two 
parameters E or F vanishes. 

Recently, on the basis of a perturbative calculation, the form of the E and F was 
conjectured to be [pOf , 

E = (o + ai)(l -B/2), F=(ao-oi)(l-S/2). (2.9) 

In fact, Corrigan and Taormina |19| gave further evidence for these formulae. They 
computed the energy spectrum of the quantized boundary states, firstly by using a 
bootstrap technique and, secondly using a WKB approximation. Requiring that the 
two methods agree with each other provides further support to the above conjecture. 
Meanwhile, they did the non-perturbative calculations in the significantly more difficult 
case when a ^ <j\. 



3 Low order perturbation theory 

Since perturbative expansion will be developed around an x-dependent static back- 
ground field configuration, it is possible to use standard Feynman rules in configuration 
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space. We will find that the theory needs three- and four-point couplings which depend 
upon the x-dependent static background. In other words, the interaction vertices in 
the Feynman diagrams will be position dependent. In order to derive the three and 
four-point couplings, we need to expand the bulk and boundary potentials in terms of 
the coupling constant (3 around the background solution. We obtain the bulk three- 
and four-point couplings: 

c£k = ^sinh(v^>o) (3.1) 



^L = ^ 2 cosh(V2/30 o ), (3.2) 



where 0o is the background solution to the equation of motion of the model. Similarly, 
the boundary three and four point couplings may be found as: 

Y (3) _ V^P ( Rfo/sft _ „ „-*/V2' 



CZndary = - (3.3) 



Cl! ndary = ^ (*ie?*°W + a e~^) . (3.4) 

Meanwhile, we may consider the linear perturbation of the equation of the motion and 
the boundary condition around the background field R IB| to yield 

, ,- 1 _|_ p 2(x-X ) 

e^o/V2 = L + e (3 5) 

I _ e 2(x-xo) ' v ' 

where Xq is related to the boundary data and satisfies 



cothxo = \ —— — -• (3.6) 
V l + o-i 

So, after some manipulation, the bulk three- and four-point couplings become 

C£L = cosh 2{x - x ) (coth 2 2{x - x ) - l) , (3.7) 



and similarly, 



^k = ^ 2 (2coth 2 2(x-x )-l) (3.8) 

^boundary = K C ° th X ~ ^0 tanhx ) , (3.9) 

y(4) _ 

48 

We also need to find the Green's function for the theory. The two-point propagator 
for the boundary sinh-Gordon model is known |16| to have the form 

+ KJ(-k,x)f(-k,x')e- ik ^ +x 'A , (3.11) 



C boundary = 7^ C ° th X + ^ taIlh X o) • (3. 10) 
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where 

/(M) = i *~ 2 ^ g ~ g ° ) (3-12) 
and represent the classical reflection factor of the model which is given by 

= / (2fc) 2 + 2ik^T+V ^T+^ 1 + 2(a + <n) \ f ik-2 \ 
c \(ik) 2 -2ikVTT^VT+^+2(a + a 1 )J \ik + 2j ' l ' j 

Indeed, the quantum reflection factor (2.8) reduces || to the classical reflection (3.13) 
in the classical limit (3 — > provided E(ao, a±, 0) = a + a\ and F(ao, o"i, 0) 

In this paper it is intended to calculate the quantum corrections to the classical 
reflection factor of the sinh-Gordon model at one loop order and up to the second 
order in the difference of the boundary parameters. Hence, it is necessary to expand 
the bulk and the boundary couplings. In fact, 

3 l + <7i 

^( 4 ) _ 1 o2 f-t , 1 ^ „4x 



"I'/'W*- 1 ' (3 - 15) 



where e = cr — a 1 . Similarly, 



^boundary— ^ \ -f ^ ^ "' (3.16) 



and 

c£L ta , = s( a " + e -4(iT^ < ') + -- (3 ' 17) 

Now let us find the expansions of /(/c, x) and the classical reflection factor K c , both of 
them appear in the Green's function (3.11). These can be shown to be 



and 



ik + 20i 2ik 
K c = - — + 



ik — 2(T\ {ik — 2a\ 



\2 



1 ik(ik 3 - Ak 2 - 6k 2 ^ - Mka x + 8af - 16 - I601) 2 
+ 2 (1 + ax)(ik - 2){ik + 2){ik - 20i) 3 C ^ ' ' 

or in a compact form 

K = K Q + K l e + K 2 e 2 + .... (3.20) 

Here, K is the classical reflection factor when the two boundary parameters are equal. 
To calculate the one loop (0(f3 2 )) quantum corrections to the classical reflection factor, 
we use the idea introduced by Kim and developed by Corrigan flB]] . In other words, 
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after perturbative calculation of the propagator the 0(/3 2 ) corrections to the classical 
reflection factor can be made by looking at the coefficient of e~ lk ^ x+x ) in the residue of 
the on-shell pole in the asymptotic region x,x' — ► — oo. 

The computations corresponding to the one loop quantum corrections to the classi- 
cal reflection factor may be performed by the standard perturbation theory which has 
been generalised to the affine Toda field theory [p^ ,[23|-p5| in the presence of an inte- 
grable boundary. There are three basic types of Feynman diagrams which contribute 
to the two-point Green's function of affine Toda field theory in one loop order. These 
are shown in figure 1. Now by looking at the three-point and four-point couplings, 
we realize that all types of these diagrams are concerned in our problem ( up to the 
second order in the difference of the boundary data). In what follows we calculate 
the contributions of types I and III Feynman diagrams to the reflection factor. The 
remaining diagrams will be treated elsewhere. 




Figure 1: Three basic Feynman diagrams in one loop order. 



4 Type III (boundary-boundary) 

Clearly, the type III Feynman diagram includes four different configurations depending 
on the fact that the interaction vertices to be located in the bulk region or at the 
boundary. This section deals with the calculations corresponding to the contribution 
of the type III Feynman diagram to the reflection factor when both vertices are placed 
at the boundary. The associated contribution may be given by 

1 ' dtdt'G{xi,t l ;0,t)G{0,t;0,t')G{0.t';0,t')G{0,t;x2,t 2 ). (4.1) 



4 (1 + a x 
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Here, the related three point couplings and the combinatorial factor have been consid- 
ered. It is evident that in this case the two-point Green function has the simplest form 

as 

G( X , t; = i / / ^:; { -ii P + *»M-^>) • (4-2) 
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where 

ik + 2(7i 



Moreover, the loop propagator has been found in Ref. flTjJ and so, 

G(0,t';0,t') = -a 1 ^^. (4.4) 
sm aiTi 

Now it is convenient to calculate the time integral of the other middle propagator in 
(4.1), i.e. 

dt'G{0,t;0,f) (4.5) 



which is equal to 



1 1 ^ (i + , . (4. 6) 

' ' 2vr 2tt cj" 2 - A;" 2 - 4 + zp \ life" - 2a x ] 




Integrating over t' gives us a Dirac delta function so, the above relation reduces to 

dk" ( 1 \ / 2ik" 



2tt \-k" 2 -AJ \ik" -2(Ji 



(4.7) 



Hence 



^(0,^0,0 = - (4.8) 
2(1 + a x ) 

Therefore, up to now the contribution (4.1) has the form 
i/? 2 aiCOsai7r e 2 f f f du dk e -M*i-*) r/ , n 



8 sinai7r (1 + cosai7r) 3 J J J 2n 2n lu 2 — k 2 — 4 + ip 

dJ dk 1 e -i"'(t-t2) 




2tt 2tt a/ 2 - A;' 2 - 4 + ip 



(e~ ik ' X2 + ^(JfeOe"**'* 2 ) • (4.9) 



First of all, it is understood to perform the transformation k — * —A; in the first term 
of the first propagator. Secondly, after integration over t which ensures energy con- 
servation at the interaction vertex, the result will be a Dirac delta function which 
immediately gives rise to the substitution u = u' . Finally, the momenta of the two 
propagator can be integrated out by taking the contours to be closed in the upper 
half-plane and considering the pole at k = k = k' = yu 2 — 4 and ignoring the other 
pole i.e. — 2%o\ (when o\ < 0) due to the fact that its residue will vanish in the limit 
Xi,x 2 — > — oo. Thus the type III (boundary-boundary) contribution becomes 

i(3 2 ai cos a l7 i e 2 f du __ Mtl _ t2) _ifc( Xl +* a ) 1 



8 sinaivr (1 + cosaivr) 3 J 2vr {ik-2a l ) 2) 
where k = 2 sinh 9. 
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5 Type III (boundary-bulk) 



Now let us calculate the contribution of the type III Feynman diagram to the reflection 
factor when the interaction vertex corresponding to the loop is placed inside the bulk 
region and the other vertex is at the boundary. As before, the combinatorial factor 
of the diagram should be considered as a coefficient factor. Moreover, in this case we 
must take into account the bulk three-point coupling in the loop vertex as well as the 
boundary three-point coupling in the other vertex. The contribution this time is 



(1 + cos ai7r) 5 



dx J J dtdt'Gix^ti^^G^^x' ,t') 

G(x',t';x',t')G(0,t;x 2 ,t 2 )e 2x '. (5.1) 



In fact, in the previous section we found the following relation which is some part of 
the contribution (5.1): 



/ 



dtG( Xl , h; 0, t)G(0, t; x 2 , t 2 ) = - I ^ e -Mti-t 2 ) e -^ 1+ * 2 ) 1 (5 . 2 ) 

J 2n (ik - 2o"i) 2 



Clearly, in the boundary-bulk contribution (5.1), it is seen that the t' dependence 
is involved only in the propagator G(0,t,x',t'). So, it is convenient to compute the 
following relation 



Jdt'G^t-x',t')=iJdt' J 



-W-*) \ ( 2ik' \ ifcV 



to' 2 - k' 2 - 4 + ip J \ik' - 2o"i 

(5.3) 

Integrating over t' generates a Dirac delta function which allows the substitution uj' = 
and hence, 

■ W y, t >) = i[%(£^)(J^). (5.4) 



Therefore, we obtain 



/ 



dt'G(0, t; x', t') = — — : -e 2x '. (5.5) 

V ' ' ' ^ 2(l + (7i) K ' 



Now let us calculate the loop propagator which has the form 

GW , , ; = i / / %% ^_J_ i + ip (l + ^.(*>— ) . (5.6) 

In fact, the above integral is divergent. However, this divergence can be removed by an 
infinite renormalization of the mass parameter in the bulk potential. In other words, 
a minimal subtraction of the divergent part can be made by adding an appropriate 
counterterm to the bulk potential. Meanwhile, if we integrate over uj" then, we will 
find 

„/,/.„/^_ 1 f dk " 1 i*" + 2<7i _ 2 ^ v 



G(x',t';x',t') = - / — - ' W1 e - 2,fc x . (5.7) 

1 ' ' ' ; 2 J 2tt ^WTZik" -2^ 1 ; 
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Hence, the remaining part of the contribution (5.1) is simplified as (apart from the 
coefficients) 



// 



dt'dx'G(0, t- x', t')G{x', t'; x', t')e 2x ' 
i 



" dx > f 1 f + 2gl W (5.8) 



4(1 + a ± ) J 2tt ^k" 2 + 4 ik" - 1o x 

In order to integrate over x', one may use the following formula 



o 



where the quantity r must be positive. The k" integration can be evaluated by closing 
the contour in the upper (lower) half-plane depending on whether o\ > (<7i < 0). 
However, due to the branch cut the contour has to run around the cut line. Note, 
the cut line stretches from k" = 2i to infinity along the imaginary axis when o\ > 0. 
Otherwise, the cut line run from k" = —1i to — oo along the imaginary axis when 
<7i < 0. So, after performing the required integrations 



// 



dt ; dx'G(0, t; x', t')G(x', t'; x', t')e 2x ' 

i ( 27raiCosai7r l + cosai7r 



167r(l + cos ain) \(1 — cosai7r) sinai7r 1 — cosai7r 



(5.10) 



Therefore, the contribution of the type III (boundary-bulk) diagram to the reflection 
factor is 



i(3e ( 27raiCosai7r l + cosai7r 



/ 



>7r(l + cos ai7r) 3 \(1 — cosai7r) sinai7r 1 — cosai7r 
2vr {ik-2a 1 ) 2 ' 

6 Type III (bulk-boundary) 

Let us consider the type III Feynman diagram when the loop vertex is at the bound- 
ary and the other interaction vertex is inside the bulk region. So, in this case the 
contribution is 

— -re 2 [ dx [ [ dtdt'G(x 1 ,t 1 ;x,t)G(x,t;0,t') 

(1 + cosaivr) 2 J ^ J J 

x G(0, t'\ 0, t')G(x, t; x 2 , t 2 )e 2x . (6.1) 

In fact, we have obtained the two middle propagators in the previous sections. 
Hence, 
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and 



/ 



. ,. cosai7r 

G O, t'; 0, if) = -ai — 6.2 

sin a\ix 

dt'G(x,t;0,t') = - * e 2x . (6.3) 
2(1 + ax) 

So, the contribution (6.1) can be shown in detail as (after doing the integration over t) 

i(3 2 e 2 dxcosaxft 
(l + cosai7r) 3 sinai7r 

x f dx [ f ^ ie~^-^ , ik(xi - x) g (fc ) e -tt(«i+«)) 
' ' ' 27T27TC,; 2 - k 2 -A + ip v ov J ' 



— CO 



x f —— _i_ f e iV(x-*2) + X (A;')e- ife '^ +3:2 ^ e 4;E . (6.4) 

J 271 ijj 1 — k' 2 — 4 + ip V / 

The integration over x may be done by using the device (5.9). Meanwhile, it is necessary 
to perform a transformation k — > —A; in the first term of the first propagator. Therefore, 
the above expression becomes 



if3 2 e 2 axcosaxn 
^l + cosai7r) 3 sinai7r 

dw dk dk' c _ iui(ti _ t2)c _ i(kxi+k , X2) 




2tt 2tt 2tt uj 2 - k 2 - 4 + ip uj 2 - k' 2 - 4 + ip 

\ k + k' -Ai k-k'-U -h + W-M -k-k 1 -Ai]' 1 ' ' 

Clearly, all that remains is to integrate over the momenta k, k'. This can be achieved by 
completing the contours in the upper half-plane and picking up the pole at k — k — k' — 
\/uj 2 — A and regarding the fact that all the other poles have no contributions because 
their residues yield exponentially damped terms in the asymptotic region xx,X2 — > 
— oo. Hence, the contribution of the type III (bulk-boundary) Feynman diagram to the 
reflection factor is 

8 2 e 2 axcosaxn f du^^-^-ik^+x,) 1 



(l + cosai7r) 3 sinai7r J 2n (2k) 2 

x {-J— + U (k) - ^—K 2 Ck)\ . (6.6) 
{2k - Ai 2 2k + Ai J 



7 Type III (bulk-bulk) 

In this section we study the type III (bulk-bulk) Feynman diagram which contribute to 
the quantum correction to the classical reflection factor. It is evident that in this case 
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both interaction vertices are located inside the bulk region. So, in this configuration 
the corresponding contribution is given by 

-16- -rre 2 / dx I dx 1 ! I dtdt'G(x 1 ,t 1 ;x,t)G(x,t;x',t') 

(1 + cosai7r) 2 J J 

G(x', t'; x', t')G(x, t; x 2 , t 2 )e 2x e 2x ' . (7.1) 

As it was shown in the type III (boundary-bulk) case, the loop propagator can be 
simplified to obtain 



ik '+ 2 °A e -*te*. (7 .2) 



4- 4 \ik\ — 2<Ji 

Moreover, the time variable t' appears only in the other middle propagator. Therefore, 
it is appropriate to calculate the following relation 



/ 



dt'G(x,t;x',t')=i I I I dt'^^ //9 6 " //9 ' ' , (e ik "^ 
111 2tt 2tt lu" 2 - k" 2 - 4 + ip 1 




+K {k")e- ik "( x+x '^ . (7.3) 



The integration over t' gives us a Dirac delta function which replaces uj" by zero and 
hence, 



/ 



dt'G(x, t- x', t>) = J ^ [-Z^TZTl) z- lk " X ' {f k " X + K,{k")e^" x ) . (7.4) 



Let us rewrite the contribution (7.1) in the expanded form in order to find how we 
may carry out the required integrations ( first setting k — > —k in the first term of the 
first propagator) 

i6 n p ? r dx r d^ fdt r i^ d ± 2 -;: ( 7 << 7 - 

(1 + cosaxTr) 2 J ^ J J J 2tt 2tt uj 2 - k 2 - 4 + ip K 

+K Q {k)e~ lk{xi+x) ) 

x f ^ J e- ik " x ' U k " x + K (k")e~* k " x ) e 2x e 2x ' 
2j2n y /kJ+Iik 1 -2<j 1 



x 




dw' dk' ie-^-^ 



2tt 2tt uj' 2 -k' 2 -4 + ip 



_ fik'{x-x 2 ) + Ko (k')e- ik ' {x+X2 A . (7.5) 



The starting point is to perform the t integration which immediately removes the 
energy variable uj'. Moreover, by multiplying the first and the third propagator then 
we obtain four pole pieces. If we perform the calculations corresponding to one of them 
( for example the first one which involves e l( - k+k ^ term), then the calculations for the 
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remaining three pole pieces may be done in the same manner, except that k + k! is 
replaced by one of k — k', —k + k! and —k — k! . So, in what follows we continue the 
computations in detail only for one pole piece. Now, the integration over x' is achieved 
by means of the device (5.9) and gives 



i 

j / (2-i(fc"+2fci))x' _ 1 ( 7 fi N 

dX€ ~k» + 2k 1 + 2i [ } 



Similarly, the x integration can be done and yields 

f dxe (2+ ^ k+k '» x (e ik " x + K (k")e- ih " x ) 

J —CO 



-K (k"). (7.7) 



k + k! + k" - 2i k + k'- k" - 2i 
So, we find that we must solve the following integral 

(3 2 e 2 f f fdudkdk' e -M*i-fc) e -i(kx 1+ k'x 2 ) 



16- 



1 + cos ai7r) 2 J J J 2ir 2ix 2tx lo 2 - k 2 - 4 + ip uj 2 - k' 2 - 4 + ip 



1 f f dk" dk x 1 /ifci + 2<Ti 
x - 



2 J J 2tt 2tt + 4 \ih - 2a t J \k" 2 + 4 J \-k" - 2k x - 2% 



X ( t + + + t + y-y-2i g ° (t,) ) ' (7 ' 8) 

It turns out to be more convenient to integrate over k" then afterwards over k\. Hence, 
first of all, we need to evaluate 



/ 



(U{ " ' ' / 1 + t ; „ * ^ o(n)- (7-9) 



2tt /c" 2 + 4 -ib" -2k l -2i\k + k' + k" -2i k + k' - k" - 2i 



The above integral may be solved by closing the contours in the upper or lower half- 
plane in order to become free of all extra poles except ±21 Therefore, after some 
manipulation we obtain 



dk" i 1 



2tt k" 2 + 4 —k" -2k l -2i\k + k' + k"-2i k + k' - k" - 2% 

1 f 1 - (7i 1 1 1 1 



Ai yl + o 1 2ki + Uk + k' — Ai 2k x + Aik + k'-2k x - Ai 

1 1 



k + k' - Ai k + k' -2ki-Ai 



(7.10) 



Next, we need to calculate the integration over k±. In order to calculate the loop 
momentum k\ integral, we will have different kinds of integrals however, all of them 
can be solved similarly. For example, consider the integral 

dkl 1 (7.11) 



y/kJ+Ak + k' -2k ± -Ai 
12 



Let us close the contour in the upper half-plane and due to the branch cut which 
extends from 2i to infinity along the imaginary axis, the contour has to turn around 
the cut line. Therefore, the integral (7.11) reduces to 

2/°°-* " (7-12) 

J 2 ^/fZT^k + k' -2iy-Ai 

and the above integral gives the following result 
dki 1 



/ 



J 1 + K fc + fc '- 4 ^W^ + 4 : (7 . 13) 



Now we can write down the solution of the k\ integration as 
1 f dh 1 fih + 2aA f 1 - a x 1 1 



Bi J 2tt + 4 \ih -1a x ) 1 1 + 01 2k x + Ai k + k' - Ai 
1 1 1 

+ 



2fci + Ai k + k' - 2k x - Ai k + k! - Ai k + k' - 2k ± - Ai 
k + k! — 2i — 2i cos ayn 



167r(l - cosaivr)(A; + k')(k + k' — Ai) 
aiCosai7r 2k + 2k' + Ai cos aire — 12i 



16 sin 3 ai7r (k + k' — Ai)(k + k' + Ai cos - Ai) 
1 k + k! — 4icosai7r — Ai 



An(k + k')(k + k' - Ai) k + k' + Ai cos a x ix - 4i ^/ (k+k'-af | ^ 

f 1 + i(A; + A;' - 4i) + f M k + k '-W + 4 ] 
x In ^ 2V . 4 = I . (7.14) 

(l + l(k + k'-A i )-^^ r ^ + A\ 

Finally, regarding the type III (bulk-bulk) contribution (7.5), all that remains is to 
integrate over the momenta k, k' . This can be achieved by means of the contours in 
the upper half-plane and considering the pole at k — k — k! — \JuP- — A and ignoring 
all the other extra poles because their contributions will vanish when X\,x<i — > —00. 
Hence, the relation (7.5) has the solution 



l + cosai7r) 2 J 2tv (2k) 



2 



k — i — i cosai7r ^aiCosai7r k + i cosai7r — 3i 



7r(l - cosai7r)/c(2/c - Ai) sin d ai7r (A; - 2i) (k + 2i cos — 2%) 
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k — 2i cosaiTT — 2i 



7rk(k — 2%) k + 2% cos — 2% ^ _ 2^2 _|_ 4 
x In 



l + i(fc-2i) + | ) /(fc- 20 2 + 4 




l + ±(fc-2i) - ^(k-2i) 2 + A 
iai cos ai7r(cos ai7T — 3) i(l + cos ai7r) 



27r(l — cosai7r) 2 2 sin 3 ai7r(l — cos aiir) ' 127r(l — cos ai7r) 
fc + i + « cosa^ aiCOsai7T — 2 cosai7r + 3i 



n(l - cos a^k (2k + Ai) sin 3 ai7r (k + 2i)(k - 2i cosaiTr + 2i) 
1 k + 2i cos ai7r + 2i 1 



nk(k + 2i) k — 2i cos ai7r + 2i 2^)2 4. 4 




,i-t(* + «) + W(* + «)'+S |: ( , 15) 

l-i(fc + 2i)-j v '{fc + 2i) 2 + 4 

Notice that the above relation depends on the In terms in a manner which is very 
inconvenient for a comparison with Ghoshal's formula. Fortunately, these terms will 
be cancelled by matching terms in the contribution corresponding to the type I (bulk) 
Feynman diagram. This is one of the interesting results which may be found in this 
paper. So, the type III (bulk-bulk) contribution takes the simple form 

P 2 ,2 f ^ U „-iu}{ti-t 2 ) p -ik{xi+x 2 ) L_ 




e / — e 



l + cosai7r) 2 J 2tt [2kf 



k — i — i cos a\Ti ^aiCOSOi7r k + % cosa^ — 3i 




7r(l - cosai7r)/c(2A; - Ai) sin 3 ai7r (k - 2i) (k + 2% cos - 2i) 



%a\ cosai7r(cosai7T — 3) i(l + cos ai7r) 



27r(l — cosai7r) 2 2 sin 3 ai7r(l — cos a^) ' 127r(l — cos ai7r) 
+ * + i cos a\ix aicosaiii k — i cos aiii + 3? 



7r(l -cosai7r)/c(2A; + 4i) sin 3 ai7r (A; + 2i) (k - 2i cos ai7r + 2i) 



(7.16) 



8 Type I (boundary) 



In fact in the case of the type I Feynman diagram, there are two possible configu- 
rations depending on the fact that whether the interaction vertex is located at the 
boundary or inside the bulk region. From now on we call these two configurations the 
type I (boundary) and the type I (bulk) respectively. Considering the boundary four 
point coupling (3.17) and the associated combinatorial factor as well, the contribution 
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corresponding to the type I (boundary) diagram is described by 
- ^ (2^ + e - 4( ^ + + ^ )2 /) j°° dt"G(x, t; 0, t" )G(0, f"; 0, t")G(0, f ; *')■ (8-1) 

Let us find the appropriate form of the two-point Green function which will be used 
many times throughout this and next sections. Now by looking at the general form of 
the propagator (3.11) and considering the expansions of the classical reflection factor 
(3.20) and the function f(k,x) i.e. (3.18) as well, the required form of the two-point 
Green function will be 



r<i / I dujdk e~^-^ jf e 2 1 4x 

Lr[X,t,X,t)-l I I 2lr27TUJ 2_ k 2_ 4 + tp \[ L+ 4(1 + 0^ (tk + 2) 6 




r 2 i 

g4x' \ ^ik(x-x') 



4(1 + a 1 f (ik-2) 



/ e 2 1 e 2 1 

+ [ K ° " 4(1 + ^ W^f K " ~ 4(1 + ^ (^2) e4l '' K ° 

+eK! + e 2 if 2 ) e-^C^') J . (8.2) 
So, in our problem the loop propagator takes the following form 

G(0 , 0, t") = / / ^^_ p jL_ (! + if. + ^ + 

e 2 1 e 2 1 e 2 1 \ 

+ 4(l + <n) 2 (ifc" + 2) ~ 4(1 + a x ) 2 (ifc" - 2) ~ 2(1 + a x ) 2 (ik" - 2) ^7 ' (8 ' 3) 

First of all it is necessary to perform a minimal subtraction in order to remove the 
divergence of the above integral. Secondly, we have solved the integral in Ref. |2(J up to 



the first order in e . Hence, what we need to do is to solve the part which is proportional 
to e 2 including four terms. Indeed, three of them can be simply manipulated and 
because of this reason we solve the first one in detail which is given by 

du " dk " 1 *K 2 {k») (8.4) 




2tt 2tt lu" - k" 2 -A + ip 



or 



1 f dk" 1 2 



2 J 2tt VA;" 2 + 4 V 4 K - 2<n + l)(fc" + 2i) 4(1 + a^^k" - 2i) 

«7i(<7? - 2<7i - 1) 1 a\-2 1 



;i + OiXof - o\ - o v + 1) + 2«7i) (1 + 0-i)(0-! - 1) (k" + 21a,) 2 

Aioi 
{k" + 2zoi) 



5) 
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Fortunately, we have already calculated these kinds of integrals in Ref. ]20] except 
the last one which is (apart from the coefficients) 

dk " 1 (8.6) 



The above integral can be converted to a complex one and regarding the branch cut, 
it reduces to 

o r dy 1 
J 2 ^T4(?/ + 2 ( t 1 ) 3 1 ' ) 

and after some manipulation we obtain 

dk" 1 3cosai7r l + 2cos 2 a!7T 



y/k" 2 + 4(ik" - 2ai) 3 8sin 4 ai7r 8sin 5 ai7r 



-a!7T. 



Now using the above formula and the required formulae in Ref. ||20|| , we may find 
the following result 




du"dk" 1 2 

e 2 / 1 + cos a\TT + cos 2 ai7r cos ai7r(2 + cos ai7r) 



(i\ix . (8.9) 



47r y sin 4 ai7r sin 5 

Finally after doing the necessary computations and simplifications in connection with 
(8.3), we obtain the loop propagator as 

G(0,W) = -a 1 ^^ + e / COSai7F 



sinai7r \27rsin a\ix 2 sin ai7r 

9 / 2 + cos 2 ai7r 3cosai7r \ 

+e 2 j—^ ^—^ti . (8.10) 

\ 47rsin ain 47rsin din J 

Let us consider the boundary contribution (8.1). The integration over t" creates a 
Dirac delta function which means we must set u/ = u. As before in previous sections, 
the remaining integral over the momenta k and k' can be achieved by closing the 
contours in the upper half-plane and taking into account the pole at k = k' = k = 
\fup- — 4. Note, all the other poles will give damped contributions in the asymptotic 
region x,x' — > — oo. So, the type I (boundary) contribution is 

i(3 2 a x + 2 2 



x 



aiCosai7r / cosa^ a x 



sinaxvr \27rsin a\ii 2 sin ayix 



,2 



2 + cos 2 ai7r 3cosai7r 
47rsin 4 ai7r 47rsin 5 ai7r 



x / duJ c -Mt-t') c -ik(x+x') 1 
2tt (2jfc)5 
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2ik 2ik 
x { + ( 

ik - 2a l (ik - 2a x ) 2 



ik(-ik 3 + 2k 2 - Uka x + Qk 2 ^ + 16 + 16<7i) 2 



2(1 + a x )(ik - 2)(ik + 2)(-ik + 2a ± f 



e 2 }. (8.11) 



9 Type I (bulk) 

Now let us evaluate the contribution of the type I (bulk) Feynman diagram to the 
classical reflection factor when the interaction vertex is placed inside the bulk region. 
It is clear that in this case we have to consider the bulk four point coupling (3.15) 
in the interaction vertex however, the related combinatorial factor is the same as the 
boundary case. The corresponding contribution may be written as 



/oo i>0 
dt" / dx"G(x, t; x", t")G(x", t"; x" , t") 
-oo J — oo 



G(x",t";x',t') {l + ^^y 4 *")- (9-1) 
As it was shown in the previous section, the loop propagator is given by 

G(x" t"-x" t") =iff — — - 1 1 + - - e 4 *" 

^ l ' t,X,tj l J J 2tt 2n uj" 2 -k" 2 -4 + tp\ + 4(l + o- 1 ) 2 (tk" + 2) 6 

fc r a a 

-e x 



4(1 + a^ 2 (ik" - 2) 

+ {*> - wk? W^f" K ° + * + ^) ^ } ' < 9 ' 2 > 

The above integral is logarithmically divergent. Nevertheless, the divergence of the 
loop integral can be removed by performing a minimal subtraction. Moreover, doing 
the integration over u" then, we obtain 



e 2 

G(x" t"-x" t") = - e 4x " 

U{X ,t ,X ,t ) 8n{1 + ai) 2 e 



+1 f dk " 1 (K 1 r ^ K 

^2 J 2tt v ^2 T 4 ( v 2(l + <Ti) 2 {ik" -2) 

+eK 1 + e 2 K 2 y- 2lk " x ". (9.3) 

Now regarding the bulk contribution (9.1), the next step is to set k — > —k in the 
first term of the first propagator. Secondly, the integration over t" leads to the fact 
that the energy variables of the first and the third propagators are equal. Hence, the 
type I (bulk) contribution takes the following form: 



/_>"// 



du dk ie-^V-^ 



W M ~ ., A ■ - 1 + 77^—^7 77T - TTTr e 



2tt 2tt uo 2 - k 2 - 4 + ip [\ 4(1 + a^ 2 (ik + 2) 



2 / j„n I I • • ii. e a x 



2 
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e 2 1 



^4x" \ ^—ik(x—x") 



4(l + <7i) 2 (ik-2) 



e 2 4x n 1 I" dk" 1 / ' e 2 1 4x tt 

87r(l + ( 7 1 ) 2e + 27 277PTll °" 2(1 + ^)2 (ifc"-2) 6 ° 

+6K 1 + 6 2 ^ 2 )e- 2jfc ' v j 

/■ rffc' ^ / / e 2 1 4." 

X ./ 2tt J 2 - k' 2 - 4 -ip\\ ■ ^ 



A(l + a 1 ) 2 (ik' + 2) 

■ ; " I ' 



' 4x' \ ik'(x"-x') 



4(1 + <7i) 2 (itf - 2) 

+£A' 1 +e 2 A' 2 )e-'-'<-" +I '>| 

*H(T^»}. <»■«> 

The integration over x" may be performed by the formula (5.9) and therefore, (9.4) 
converts to 

-2j3 2 [ [ [ — —— e -Mt-t') e -i(kx+k'x') i i 



27T 2tt 2tt uj 2 - k 2 - 4 + ip u 2 - k' 2 - 4 + ip 

1 



[dJf__f_[ K (h „, ( 1 1 
7 27r v ^T4\ 01 ^4(1 + ^)^ 



4(1 + ax) 2 ik + 2k + k' - 2k" - 4i 
1 1 1 



4(1 + <n) 2 ik' + 2 k + k! - 2k" - Ai 



| 1 1 Mk') | 1 1 K (k) 



4(1 + ai) 2 ik + 2 k - k' - 2k" - Ai 4(1 + a x ) 2 ik' + 2k' — k — 2k" - Ai 
1 1 K (kf) 1 1 K (k) 



4(1 + a^ 2 ik' — 2k — k' — 2k" - 4i 4(1 + a^ 2 ik -2 k' - k- 2k" - 4i 
1 1 K (k)K (k') 1 1 K (k)K (k') 

' a / -1 \ Q • 7 / r\ 1 . 1 1 r\ 1 11 A ' 



4(1 + (Ti) 2 ik' -2k + k' + 2k" + Ai 4(1 + a^ 2 ik - 2 k' + k + 2k" + Ai 
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2(1 + en) 2 ik" -2\k + k' - 2k" - 4i k — k' — 2k" - Ai' 

+ 2(i + CTl )2 A ' o(r) (i + if - 2*» - 5 + r^¥^W^Ti Ko{k ' ) 

Meanwhile, there is another term which is the remaining part of the contribution i.e.: 

(3 2 e 2 f f f dw dk dk' c - <h> (f-t>) c - i(fc3!+fc v) 1 1 



2tt(1 + (Ji) 2 7 J 7 27r27r27r uj 2 - k 2 - A + ipu 2 - k' 2 - A + ip 

{k + k'-Az + k-k'-A^^ 

+ V^~u K ^ + ^iT^XoWoW) • (9-6) 

Next, we need to integrate over k" however, the corresponding calculations are 
tedious. We will discover two important facts during the computations. Firstly, all the 
terms involving the following kind of integral 

dk" 1 



Vk" 2 + Ak + k'- 2k" 

which is proportional to 9 (after doing the integrations over k and k! and using the 
fact that k = k = k' = 2 sinh 9) will be cancelled. This is consistent for a comparison 
with Ghoshal's formula. Secondly, there are many other terms which contain this type 
of integral 

dk" 1 



y/k" 2 + A k + k'- 2k" - Ai 

and all of them (after assembling) will be canceled with the counterpart terms in the 
type III (bulk-bulk) Feynman diagram as we mentioned in section seven. This result 
is also compatible with Ghoshal's formula. Meanwhile, further simplification can be 
made by using the values of k and k' given by their poles and simplifying the integrand. 
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So the type I (bulk) Feynman diagram has the following contribution 

iB 2 e 2 f iuj (t-t') -ik(x+x')_ 1 



n J 2n (2k) 2 

cosai7r diix ( 1 1 1 Kl(k) 



1 + cos aiTr) 2 sin aivr \ik + 2 ik-2 cos a x -K + 2 ik-2 ik + 2 cos g^tt - 2 

, 2ai7rcosa 1 7r\ 1 1 f 

+ 1 — 7Z— — T^— K oik) 



smai7r / sin ai7r A; 2 + 4 

. 1 ai7r 2i/c 7r . . 

+ -t + — * + K 2 (k) 

4 sniai7riA;_2cosai7r V^ + V 

1 diix cosai7r\ 2iA; r ^ /f . a\K ik r ^->,ts 

-ATi(fc) + — T? K 2 (k) 



2 sin ai7r 2 sin d ) ik - 2 cosaivr sinai7r ik + 2 cosaivr 

ai7r cosai7r(2 cosai7r — 3) / 1 K (k) -^o(^) 



2(1 + cosaiTr) sin 3 ai7r \^iA: - 2 cosa^ + 2 l-cosa^ iA; + 2 cosa^ - 2 
7 cosai7T \ 2cosai7r 

+ ( T77— 2- _ + :r-i ) K (k) - — 



12 sin ax7r 2 sin a X T\ ) ' sin a x i\ (ik - 2 cosa^ f 

^1 ^ aiix ^ cosai7r 4 cos 2 ai7r — A; 2 



4 VsinaxTr / sin aivr (jfc - 2 cosaivr) 2 

--^^(K 2 (k) + K 2 (k)K 2 (-k)) 
2smai7T V / 

1 / ai7T (1 + 3cos 2 ai7r)\ A; 2 

+ z 2cosai7T- 1 \ , — )— 9.7 

2sin 4 ai7r V 2sma 1 7r J (ik - 2 cos a^ f 

10 Ghoshal's formula up to the second order 

Ghoshal's formula (2.8) for the sinh-Gordon quantum reflection factor up to one loop 
order is given by: 



K q (9) ~ K c (9) (l - ^sinhtf F(9)j , (10.1) 



where 

m = 



cosh 9 + 1 cosh 

ei ei 



cosh6> + sin(e 7r/2) cosh6> — sin(e 7r/2) 

j ii ii no 2) 

cosh0 + sin(/ O 7r/2) coshfl - sin(/ 7r/2) ' v ' ; 

In calculating (10.2) we have made use of the expansions of E and F to 0(f3 2 ): 

B 2 B 2 
E~e + ei ^- F~/o + A^, (10.3) 
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where 



eo = clq + 0,1 



and 



fo — Oq — 0\. 



Since K q = K c + 5K C , we find that 



id 2 



Hence, expanding to 0(e 2 ), we obtain 

m ' 1 ! 



+ 



ei 



ei 



cosh + 1 cosh cosh 9 + sin ai7r cosh 9 — sin ai7r 



eiecosai7r 



1 



+ 



2sinai7r \ (cosh 9 + sin aiir) 2 (cosh 9 — sin ai7r) 2 



1 



1 



SsinaxTT \ (cosh 9 + sin g^tt) 2 

1 



eie 2 cos 2 ai7r 



(cosh 9 — sin g^tt) 2 
1 



4 sin ai7T \ (cosh 6 1 + sin ai7r) 



(cosh^ — sinai7r) 



+ 



sin ai7r cosh 6 1 



where [20 



ei 



-2ai + 



+ 0(e 2 



10.4) 



(10.5) 



10.6) 



(10.7) 



7rsinai7T 

and /i is proportional to e. 

On the other hand, the relative correction to the classical reflection factor K c is 
given by (using (3.20)) 



SK C 



■■ K^SKo + e (^Kq^ 6 K\ - K^ 2 K X 5K ) 
+ e 2 [K^8K 2 - K^ 2 K x bK x - K 2 K 2 6K + K^K 2 5K ) . (10.8) 

In fact, the corrections to K (SK Q ) and Ki(SKx) have been calculated in Ref. JT^| and 
Ref. 1 2"0fl , respectively. Clearly, all the contributions that have been obtained in this 
paper are involved in the correction to K 2 i.e. 5K 2 . It is seen that the functional form 
of the second order quantum corrections to the classical reflection factor corresponding 
to the type I and type III Feynman diagrams are consistent with Ghoshal's formula. 



11 Discussion 

In this paper we studied the sinh-Gordon model restricted to a half-line by boundary 
conditions which are compatible with integrability. We calculated the second order 
quantum corrections to the classical reflection factor of the model. The principal 
purpose of the perturbative calculation is to test the quantum reflection factor of the 
sinh-Gordon particle given by Ghoshal up to the second order in the difference of the 
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two boundary parameters. In fact, Ghoshal found the general form of the quantum 
reflection factor of the sinh-Gordon model. However, apart from two special cases 
(Neumann and Dirichlet boundary conditions) Ghoshal's formula could not provide 
a complete relationship between the reflection factor and the boundary data. As we 
mentioned before, on the basis of a perturbative calculation, the relationship between 
the various parameters was conjectured to be pO 



E=(a + a 1 ){l-B/2), F = (a - ai)(l - B/2). (11.1) 



Corrigan et.al RBI tried to give further support for these formulae. They used two 
different technique i.e. the boundary bootstrap method and the WKB semi-classical 
approximation approach. Their non-perturbative calculations have been performed for 
the complicated case of the general boundary conditions in which the bulk symmetry 
is violated. Alternatively, the conjecture (11.1) may be checked on the basis of pertur- 
bative calculations up to the second order. It is understood that this paper provides 
most part of the collection of ingredients one need. 

The computations associated with the type II Feynman diagram have not been 
performed in this article. In fact, in this diagram the middle momenta will be related 
to each other in an intricate manner due to the double Green functions. Nevertheless, 
it is expected that the same procedure will be applicable to the remaining diagram as 
well. Moreover, many formulae of this paper will be helpful. 

All the affine Toda field theories whose corresponding Lie algebras are simply laced 
or non-simply laced algebras have known exact quantum S-matrices. However, in the 
presence of a boundary the boundary S-matrices or the reflection factors of the theories 
are largely unknown. In other words, the analogue of Ghoshal's formula has not been 
formulated for most of the models. Besides, classifying the quantum integrability of 
the affine Toda field theories has not been completed. 
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